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This paper describes the application of the rain attenuation theory of 
Ryde and Rijde to the design of Radio Systems. It shows that an upper 
bound on the outage time due to rain attenuation can be computed from a 
measured point rain rate distribution. The paper also describes a suitable 
rain gauge. 

I. INTRODUCTION 

Heavy rainfall on a radio path absorbs and scatters power transmitted 
at frequencies above 10 GHz and causes large fading of received signals. 
At 20 GHz, for example, the attenuation due to a uniform rain rate of 
100 mm/hr is about 10 dB/km. Rain attenuation is so severe at these 
frequencies that for some applications transmission paths must be re- 
stricted to a few kilometers or less rather than the tens of kilometers 
common at lower frequencies. Since the cost of a radio system increases 
with the number of repeaters it is important to use the longest path 
allowed by the transmission objectives. This path length can be deter- 
mined accurately only if the fading outage due to rain attenuation can 
be predicted. 

Bussey estimated fading statistics on a microwave path from point 
rain rate data. 1 He used the rain attenuation theory of Ryde and 
Ryde 2-4 to convert rain statistics to fading statistics and since 1950 his 
results have been used in the design of radio systems. 5 However, as 
operating frequencies increase and path lengths get shorter, increased 
precision of fading estimates is required for optimum radio system 
design. 

Over the years a number of experiments have been performed in 
which attenuation was measured on a path at specific times and com- 
pared with values computed from rain rates measured by rain gauges 
spaced along the path near ground level. Here too, the theory of Ryde 
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and Ryde was used. In general there is wide disagreement between com- 
puted and measured values. Also, Medhurst questions the validity of 
the application of the theory to a practical rainfall situation. 6 Since 
the theory was derived for uniform rain the question is a good one — 
rain on a microwave path is seldom uniform. 

In this paper the theory of attenuation by uniform rain is applied to 
the practical rainfall situation. The radio path is denned as the volume 
of the first Fresnel zone and the rain attenuation is assumed proportional 
to the number of raindrops in this path volume. Of course, this expres- 
tion reduces to that of Ryde and Ryde when the rain is uniform. 

Rain rate is a vector which can be written as the product of a rain 
density and the velocity of raindrops. Rain density is proportional to 
the number of raindrops per unit volume. 

Since rain rate is a vector the expression for attenuation in terms of 
rain density in the path volume can be transformed by the divergence 
theorem into an expression for attenuation as a function of the rain rate 
on the surface of the path volume. From this formulation the following 
results emerge: 

(i) A natural definition of rain rate which is appropriate to the radio 
situation. 

(ii) A time interval, T , exists during which no significant fade can 
occur. T is determined by the path length, the frequency of operation 
and the speed of raindrops. 

{Hi) A rain gauge is described which is suitable for measuring rain 
rate in accordance with the definition mentioned in (i). 

By applying these results, an upper bound on the outage time due to 
rain attenuation is derived. The bound can be computed from a meas- 
ured point rain rate distribution using the results of uniform rain 
theory. The bound can be made tight by the proper choice of rain rate 
integration time interval. A method of estimating this interval from 
measured path loss distributions is given. 

II. GENERAL CONSIDERATIONS 

2.1 The Radio Path 

The radio link consists of two narrow-beam antennas pointing directly 
at each other over a distance of a few hundred to a few thousand meters. 
The space, or volume, of the path is taken to be the first Fresnel zone. 7 
This means that only the energy confined to that volume contributes 
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significantly to the total energy collected by the receiving antenna. 
The first Fresnel zone is a long, thin, prolate ellipsoid of revolution. 
For a path of length L at wavelength X, it has a major axis L and equal 
minor axes (XL)* and is terminated at the ends by the antennas. The 
radio path is defined as the volume enclosed by the first Fresnel zone 
and the two antennas. Figure 1 is a sketch of the path. When we speak 
of rain falling on the path we mean rain falling through this volume. 

2.2 Rain Rate as a Vector 

A theory of rain attenuation has been formulated by Ryde and 
Ryde, 2-4 and others, and a good account of it is given by Medhurst. 8 
The attenuation in a radio path depends upon the number and size of 
the raindrops and not explicitly upon their speed or direction. But the 
quantity usually measured is rain rate and it does depend on the speed 
and direction of the raindrops. Since rain rate is the product of a density 
and a velocity, it can be interpreted as a vector. 

Let there be a uniform distribution of N D drops of water per cubic 
centimeter in the space between two antennas. The drops are spherical 
with diameter D and velocity v D . The fraction of volume occupied by 
water is defined as the rain densitv 



Po - IN D D\ 



(1) 



Rain density is a dimensionless, real, nonnegative quantity. The rain 
rate for drops with diameter D and velocity v D is 

Rd — PdVd • 

The direction of the rain rate is the direction of travel of the drops. 
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Fig. 1 — Radio path. 
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The vector expression for rain rate is, therefore, 

R D = PdVd , (2) 

where a boldface letter denotes a vector quantity. 

In general a rain storm has drops of many diameters and the total 
rain rate is a summation over the drop diameters present. 

R = £ PD v D . (3) 

D 

2.3 Attenuation and Rain Density 

For a rain consisting of uniformly distributed drops with diameter D 
the attenuation of radio waves with wavelength A is 

Attenuation = 4.343 ^r~ A D X 10 5 dB/km (4) 

where A D is a function of the drop diameter, the wavelength, and the 
dielectric constant of water. Substituting from (1) the attenuation is 

a D = k(\, D)L Pl , dB, (5) 

where 

k(X, D) = 3 X 4.343 ^4$ X 10 5 . 

TV 1) 

The result in (5) is extended to the case of nonuniform spatial distri- 
bution of raindrops by replacing the uniform rain density in (5) with 
the average rain density in the radio path. The path attenuation is 
therefore assumed to be 

a D (t) = k(\, D) | fff p D (x, y, z; t) dV dB. (6) 

V 

This expression reduces to (5) for uniform rain density. 

There are neither sinks nor sources of rain in the radio path so, for 
constant drop diameter D, the hydrodynamic equation of continuity 
which applies is 8 

V-( Po v fi )+^= o, 

where the time and space variables have been omitted for convenience. 
Substitution into (6) gives 
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^P = fc(X,i»|///[_V(p^)]d7 

and transforming to a surface integral by the divergence theorem we get 
*«Jj& = k(KD)^f s [-p D v D -d6] i (7) 

where the integral is over the surface enclosing the volume V. 

Expression (7) relates point rain rate to path attenuation. The vector 
differential d6 has a magnitude equal to the differential area of the 
surface S, is normal to it and directed outward. Consider a small area 
on the surface S. The rain rate at this point is p D v D ; the component of 
this rain rate which is normal to the surface and directed into the path 
volume is — p D v n -rfd/| dd |. The integral over the surface S is the in- 
crease in water density on the path in unit time which causes an increase 
in attenuation. 

Let the surface £ enclosing the radio path volume V be described by 
orthogonal parametric curves on the surface. 10 A point P(x, y, z) on the 
surface can be written P(u, v) where, 

(XL) 1 . 

x = — x — sin u cos V 

L 

y = — cos u 



(XL) 5 . 
z = sin u sin v. 

With this transformation, and substituting R D from (2), (7) becomes 

^P = fe(X, D) \ T jf s [-Ro(u, v; t)]-ds, (7a) 

where r/s is the transformed vector surface differential. Integrating both 
sides over time T and dividing by T results in the following equation. 

(8) 

where the unit vector r is defined by R D (w, v; t) = R D (u, v; t)r. 

The left side of (8) is an approximation to the rate of change of 
attenuation since, by definition, 
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da D (t) _ .. [a D (t + T) ~ «„(fl] 
dt ~ ™ T 

Whether attenuation or rain rate is measured, the measuring instrument 
has some time constant, or integration time, T. Since, in a microwave 
system, an outage of a millisecond may be significant, an important 
question arises — can we be sure that the measurement will give us all 
the important data? In other words does an integration time T exist 
such that for T ^ T„ , (8) is a good approximation to (7a)? The ex- 
istence of T is justified by the following physical argument. 

For t « t let there be no rain on the path. At t = t let the rain rate 
at every point on the surface S be R„ = p Q v Q and be directed inward in 
the direction of the shortest line from the surface to the path axis. In 
order that a substantial fade occur, the rain will have to travel a sub- 
stantial fraction of the shortest distance from the path surface to the 
path axis. The average of this distance is (XL)*/3 and the average time 
required for the rain to reach the axis is (\L)*/3tfo • An integration time 
T « {\L)*/3v is therefore sufficient since substantial fades will not 
occur in times less than this. If T is chosen small enough, say T = T , 
(8) is a good approximation to (7 a), and can be written 



^«^^?//.{-fer" w -^^4- di - 



(9) 



The quantity in brackets is the rain rate at the point {u, v) averaged 
over time T and defines a suitable rain rate measurement. A rain gauge 
which measures rain rate in accordance with this definition is described 
in the Appendix. If the rain rate is known at every point on the surface 
of the radio path the time derivative of attenuation can be computed 
from (9). 

Since rain rate cannot be measured at all points on the surface of the 
path we consider what can be done with a more reasonable experiment — 
one or more rain gauges near the ground in the vicinity of the path. 
No satisfactory theoretical expression has been derived for computing 
attenuation from rain rate measurements made near the ground in the 
vicinity of the path. And the wide disagreement between computed 
and measured attenuations in experiments of this type described in the 
literature support the conclusion that the empirical expressions used 
are also unsatisfactory. 6 Also, the requirements of the sampling theorem 
must be met if the attenuation is to be computed from sampled rain 
rates. 11 Visual observation of rainfall reveals a spatial structure so fine 
that an unreasonably close spacing of rain gauges would be required to 
meet these requirements. 
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III. POINT RAIN RATE AND PATH ATTENUATION DISTRIBUTIONS 

3.1 Important Assumptions 

Fortunately, the attenuation as a function of time is not required 
for the design of radio links; for radio link design the fraction of time 
that the path attenuation exceeds the fading margin is the important 
parameter. Thus, we need only a suitable statistic of path attenuation 
which can be related to a similar statistic of rain rate at ground level 
in the vicinity of the path. In this section a point rain rate distribution 
function is denned and related to the path attenuation distribution 
function. 

Let the rain rate be R(w, v; t) on the surface S of the path. At some 
point near the path — on the ground beneath the path, for example — a 
rain gauge measures a sequence of rain rates given by 

1 i»i.+ (*+l)T 

R n (T, k, x, y, z) = 7n [-R(z, V, z\ t) -k] dt (10) 

where k is a unit vector normal to the collecting surface of the rain 
gauge and x, y, z are the space coordinates of the rain gauge. Suppose 
that rain rate measurements have been made for a very long time. The 
data available is a large number of rain rates R„(T, k, x, y, z), one for 
each interval T. The data is organized by choosing a rain rate R and 
computing the fraction of intervals T for which the rain gauge recorded 
rates less than R . This fraction is denoted by 

P[R n (T, k, x, y, z) ^ R„], (11) 

and is a point rain rate distribution function; it is a function of R , the 
integration time T, the location of the rain gauge, and the pointing di- 
rection k. The optimum direction for k is expected to be vertically up- 
ward in many regions; in any case, the rain gauge must be pointed in the 
direction k such that for the high rain rates of interest, that is, for 
R > R m ) 

P[R n (T, k, x, y, z) S R ] ^ P[R n (T, 1, x, y, z) < R ], 

where 1 is any unit vector. 

For the radio systems of interest, the rain rate, R m , may be chosen 
at least one order of magnitude greater than the mean rain rate. In 
this country the mean rain rate is on the order of 0.1 mm/hr whereas a 
reasonable value of R m may be 10 mm/hr. 

Let the radio path be divided into a large number, n, of volume 
elements such that the rain rate is uniform in each element. The average 
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rain rate on the path is 

R.M -ijb^fo >v< ,ii ;*)• (12) 

n ,=i 
In integral form this is written 

RUO = \ Iff R(x, y, z; dV. (12a) 

Two assumptions are made: 

(i) In a region containing the radio path the point rain rate dis- 
tribution function is independent of position and (11) can be written 

P[R n (T, k, x, y, z) S R.] = P[R n (T) ^ R ]. (13) 

(ii) For the rain rates of interest, that is, for R > R m , and for 
integration time T, the distribution function of the average rain rate 
on the path is greater than, or equal to, the point rain rate distribution 
function. 

P[R„.(t, T) ^ R m ] ^ P[R n (T) ^ R ], (14) 



where, 



RUt, T )^\\ t 5.«(0 dL 



The first assumption is that the point rain rate distribution function 
is the same whether it is measured below the path, on the path, or 
near the path. It does not mean that the rain rate at any time t is the 
same everywhere in the region — an essential distinction. The assumption 
does not mean that rain rate is a stationary random process in either 
the wide sense or the strict sense. In statistical language it means that 
the rain rate can be considered a first order stationary process over a 
small interval. 12 The second assumption reflects Bussey's observation 
that high rain rates extend over smaller areas than do lower rain rates. 

3.2 A Bound on the Path Attenuation Distribution 

If the speed of the raindrops does not change while in the radio path, 
the attenuation can be written in terms of rain rate. From (6), 

a n (t) = ^^ | /// Rn(x, y,z;t) dV . (15) 

v 

Let the rain have the Laws and Parsons distribution of drop diam- 
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eters. Then 

R D (x, y, z; t) = R(x, y, z; t)p D , (16) 

where p D is the fraction of water in the rain consisting of drops of 
diameter D. Substituting into (15) and using the definition (12a) the 
expression for attenuation is 

«o(Q = ^^LpaRUt). (17) 

The total attenuation is 

a(t) = LR avt (t) Z^T^-Po. (18) 

D V D 

The quantity represented by the summation has been computed by 
Ryde and Ryde and by Medhurst for the Laws and Parsons drop 
diameter distribution and for the terminal velocities of water drops in 
still air. 8 

In Section 2.3 we showed that negligible changes occur in a(t) in a 
time T ^ T . Thus, (18) can be written 



a(t) « a(t, T.) = LRUt, T.) E ^r^ 2 Pi, , (19) 



k(X, D) 

aw ~ a\t, i B ) = xjix av .\i, l ) £_, 

u 

where 

«(t, T) = fl* «(0 dt. 
The path attenuation for a uniform rain rate R„ is, from (18), 

a^LR.Z^^Pn ■ (20) 

n v D 

The desired bound can be found by substituting from (19) and (20) 
into (14). 

PW) ^ a )) ^ P[R n (T a ) ^ R„]. (21) 

This bound says that if the measured point rain rate distribution is 
converted to an attenuation distribution by (20), the outage time 
predicted is greater than, or equal to, the outage time that occurs on 
the path. The application of this bound is illustrated in Section V. 

IV. EXPERIMENTAL DETERMINATION OF INTEGRATION TIME 

It was shown in Section 2.3 that if the rain rate integration time T 
is short enough no significant fades will be missed; specifically if T = 
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T « (\L) h /Sv , it is small enough. This determination of T is based 
on considerations as to what could happen on a path. Fades may occur 
slowly compared with T so it is worthwhile to determine whether a 
larger integration time may be practical. 

Suppose that path attenuation and point rain rates have been meas- 
ured with an integration time T ^ T . Then the distributions com- 
puted from the measurements will satisfy the inequality (21) which 
can be written 

P[«(t, T) ^ a e ] ^ P[R n (T) g R.], T ^ T . (22) 

Now, from (14), this inequality holds for any T and since T ^ T all 
fading of significance is included. 

Let the attenuation distribution be computed from the path loss 
measurements for integration times of T, 2T, 37\ • • • , as long as the 
distribution remains substantially unchanged. The point rain rate 
distributions computed for the same integration times are such that 
the inequality holds and 

P[a(t, mT) ^ a ] ^ P[R„(mT) ^ R ], (23) 

for all m = 1, 2, 3, • • • . Now suppose that P[a(t, mT) ^ a ] remains 
substantially unchanged for all integers m up to M. Then all of the 
corresponding rain rate distributions result in valid upper bounds on 
outage time as shown by (23); one of these will be the least upper 

bound. 

From experiments of this kind, practical values of integration time 
can be determined. It is expected (but not proven) that the integration 
time which results in the least upper bound will be the largest value 
for which the attenuation distribution remains unchanged, that is, MT. 

V. DISCUSSION 

Experimental rain rate distributions which meet the requirements 
of this theory are not available. For this reason a careful experimental 
verification cannot be made now. There is reason for optimism, however. 
The upper bound on outage time computed from Bussey's 1 one-minute 
point rain rate distribution is remarkably close to the one-minute 
attenuation distributions reported by Semplak and Turrin. 13,14 These 
distributions are shown in Fig. 2. Semplak and Turrin also report that 
the attenuation distribution remains unchanged for shorter integration 

times. 14 

The simplest application of this theory to the design of a radio path 
requires only a point rain rate distribution measured as described in 
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Fig. 2 — Comparison of computed and measured attenuation distributions for 
one minute integral ion time. 



Section III. In this case the rain rate integration time is computed 
from the frequency, the path length, and the maximum raindrop veloc- 
ity. This measured point rain rate distribution can be converted to 
an upper bound on attenuation outage time by use of expressions 
(20) and (21). The upper bound obtained this way is not necessarily 
the least upper bound whereas the least upper bound is required for 
optimum radio system design. The least upper bound can be computed 
from the point rain rate distribution if the corresponding optimum rain 
rate integration time is known; the optimum rain rate integration time 
can be determined from a path attenuation experiment as described 
in Section IV. 

It is important, then, to determine the optimum integration times 
in those regions of the country where radio systems above 10 GHz may 
be used. The optimum integration time is a function of wavelength, 
path length, and the climate in which the path is located. A few path 
loss measurements, located in different climatic regions, would yield 
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valuable results. In these experiments attenuation distributions would 
be measured as functions of path length and wavelength, and optimum 
integration times computed from this data and the measured point rain 
rate distributions. It may be, for instance, that the optimum integration 
times are about the same everywhere; if so, the path loss experiments 
would show it and, thereafter, only point rain rate measurements 
would be required. 

The accuracy with which the least upper bound predicts the outage 
time on a radio path cannot be stated precisely until further experi- 
mental data is available. It may be anticipated, however, that the 
accuracy will decrease as the path length increases. For example, if 
the path is long compared with the dimensions of thunderstorms, the 
least upper bound prediction will probably be pessimistic, especially 
for large fades. On the other hand, for paths shorter than the dimensions 
of thunderstorms the least upper bound prediction may be accurate. 

A rain gauge, from which the rain rate in each fixed integration 
interval can be determined, is suitable for the determination of the 
point rain rate distribution function defined in Section III. Morgan 
has built and described such a rain gauge recently, and another rain 
gauge proposed for this purpose is described in the Appendix. 

No attempt has been made to include the effects, on attenuation, of 
raindrop distortion, temperature, radio wave polarization, and so on, 
in this theory. 
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APPENDIX 

An Instrument for Measuring Point Rain Rates 

The instrument described measures rain rates in accordance with 
the definition of (9). Figure 3 illustrates the basic element in the rain 
rate instrument — a depth gauge consisting of a funnel, a cylindrical 
capacitor, and a shutter for draining the capacitor. The funnel is exposed 
to the rain for a specified time T and then covered. The rain which 
passed through the collecting area of the funnel drains into the cylin- 
drical capacitor and forms a column of water as shown. The dielectric 
constant of water increases the capacitance — the greater the height of 
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Fig. 3 — Capacitor depth gauge. 



the water column, the greater the capacitance. When the funnel has 
completely drained, the capacitance (and hence the volume of water) 
is measured. 

There are two features of importance about this depth gauge. 

(i) The interval of exposure, or integration time, can be determined 
precisely by careful operation of the shield over the funnel, and is not 
dependent on the rain rate. 

(ii) The time allowed to drain the funnel and measure the capacitance 
is independent of the exposure interval T. Sufficient time can be allowed 
to drain the funnel and stabilize the water column prior to measuring 



134 



THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1970 



the capacitor. This eliminates fluctuations due to the random behavior 
of water flow on the surface of the funnel. 

When the shield is over the depth gauge no rain is collected. The 
complete rain gauge consists of a number of depth gauges arranged as 
shown in Fig. 4. In this illustration, ten depth gauges are shown, with 
number 2 in position to collect rain. When the interval ends, the shield 
rotates, covering depth gauge number 2 and exposing number 3. 

A rotating shutter for draining the depth gauge capacitors is shown 
in Fig. 4B. The shutter is fixed to a common shaft with the funnel 
shield of Fig. 4A. The operating sequence is as follows. There are ten 
time intervals of length T in a single rotation of the shaft. With the 
aid of Fig. 3 the following sequence can be seen to occur. 



Time Interval 
Number 
1 
■2 
3 
4 
5 
(5 
7 
8 
9 
10 



Status of Depth Gauge 



Collecting Water 
Draining Funnel 



Measure Capacitance 
Draining Capacitor 
Draining Capacitor 



#:5 



Draining Capacitor 
Collecting Water 
Draining Funnel 



Measure Capacitance 
Draining Capacitor 



The rain shield steps rapidly from gauge to gauge. There is always 
one gauge collecting rain; one measurement is made in each time in- 
terval T. If the measurement starts at time t and the rain gauge points 
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-ROTATING 
SHIELD 




ROTATING SHUTTER 
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NO. 3 AND 4 DRAINING 



CAPACITOR 

CONNECTIONS 

TO 

MEASURING 

CIRCUIT 



(b) 



Fig. 4(a) Top view of rainrate gauge, (b) Section of rainrate gauge. 
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vertically upward, the output of the instrument is a sequence of meas- 
urements 

«.«, *, y, *) = ip / R(t) dt. 

REFERENCES 

1. Bussey, H. E„ "Microwave Attenuation Statistics Estimated from Rainfall 

and Water Vapor Statistics." Proc. I.R.E., 3S, No. 7 (July 1950), pp. 781- 
785. 

2. Ryde, J. W., "Echo Intensity and Attenuation due to Clouds, Rain, Hail, 

Sand and Dust Storms at Centimeter Wavelengths," Rep. 7831, General 
Electric Company Research Laboratories, Wembley, England, October 
1941. 

3. Ryde, J. W., and Ryde, D., "Attenuation of Centimeter Waves by Rain, 

Hail and Clouds," Rep. 8516, General Electric Company Research Lab- 
oratories, Wembley, England, August 1944. 

4. Ryde, J. W, and Ryde, D., "Attenuation of Centimeter and Millimeter 

Waves by Rain, Hail, Fogs and Clouds." Rep. 8670, General Electric Com- 
pany Research Laboratories, Wembley, England, Mav 1945. 

5. Hathaway, S. D.. and Evans, H. W., "Radio Attenuation at 11 KMC," 

B.S.T.J., 8S, No. 1 (January 1959), pp. 73-97. 

6. Medhurst, R. G., "Rainfall Attenuation of Centimeter Waves: Comparison 

of Theory and Measurement," IEEE Trans. Antennas and Propagation, 
AP-13, No. 4 (July 1965). pp. 550-564. 

7. Slater, J. C, and Frank, N. H, Introduction to Theoretical Physics, New 

York: McGraw-Hill, 1933, pp. 307-314. 

8. Petterssen, S.. Weather Analysis and Forecasting, New York: McGraw-Hill, 

1956, pp. 7-8. 

9. Lass, H., Vector and Tensor Ajialysis, New York: McGraw-Hill, 1950, p. 114. 

10. Struik, D. J., Differential Geometry, Cambridge, Massachusetts: Addison- 

Wesley Press, 1950, p. 64. 

11. Ruthroff, C. L., "Microwave Attenuation and Rain Gauge Measurements," 

Proc. IEEE. 57, No. 6 (June 1969). pp. 1235-1236. 

12. Papoulis, A., Probability, Random Variables, and Stochastic Processes, New 

York: McGraw-Hill. 1965. 

13. Scmplak. R. A., and Turrin, R. H, "Some Measurements of Attenuation bv 

Rainfall at 18.5 GHz," B.S.T.J.. 4S, No. 6 (Julv-August 1969), pp. 174S- 
1756. 

14. Semplak, R. A., unpublished work. 

15. Morgan. B. E.. "A Rainfall Rate Sensor." Technical Memorandum ERLTM- 

NSSL 42, Experimental Science Sen-ices Administration Research Lab- 
oratories. National Severe Storms Laboratory, Norman, Oklahoma, No- 
vember 1968. 



